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ABSTRACT 

We describe the ROSAT All-Sky Survey — Center for Astrophysics Loose Systems (RASSCALS), the 
largest X-ray and optical survey of low mass galaxy groups to date. We draw 260 groups from the 
combined Center for Astrophysics and Southern Sky Redshift Surveys, covering one quarter of the sky 
to a limiting Zwicky magnitude of m z — 15.5. We detect 61 groups (23%) as extended X-ray sources. 

The statistical completeness of the sample allows us to make the first measurement of the X-ray 
selection function of groups, along with a clean determination of their fundamental scaling laws. We 
find robust evidence of similarity breaking in the relationship between the X-ray luminosity and velocity 
dispersion. Groups with a v < 340 km s _1 are overluminous by several orders of magnitude compared to 
the familiar Lx oc c 4 law for higher velocity dispersion systems. An understanding of this break depends 
on the detailed structure of groups with small velocity dispersions a v < 150 km s . 

After accounting for selection effects, we conclude that only 40% of the optical groups are extended 
X-ray sources. The remaining 60% are either accidental superpositions, or systems devoid of X-ray 
emitting gas. Combining our results with group statistics from N-body simulations, we find that the 
fraction of real, bound systems in our objectively selected optical catalog is between 40%-80%. 

The X-ray detections have a median membership of 9 galaxies, a median recession velocity cz = 7250 
km s _1 , a median projected velocity dispersion a v = 400 km s" 1 , and a median X-ray luminosity 
L x = 3 x 10 42 

^ioo er § s 1 ' where the Hubble constant is H n — 100/iion km s 1 Mpc x . We include a 
catalog of these properties, or the appropriate upper limits, for all 260 groups. 
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1. INTRODUCTION 

The ROSAT and ASCA missions have shown that even 
low mass systems of galaxies contain a hot intergalac- 
tic plasma. Many of the Hickson (1982) compact groups 
(HCGs) are embedded in diffuse X-ray emission detectable 
by the ROSAT (Ebeling, Voges, & Bohringer 1994; Pildis, 
Brcgman, & Evrard 1995; Ponman et al. 1996). Other 
ROSAT studies of smaller group samples (e.g. Henry et 
al. 1995; Burns et al. 1996; Mulchaey et al. 1996; Mah- 
davi et al. 1997; Mulchaey & Zabludoff 1998) confirm the 
existence of an intergalactic plasma with an average tem- 
perature kT w 1 keV in many loose groups. 

Although heterogeneous studies abound, an objective 
survey of the nearby universe for X-ray emitting groups is 
lacking. A search based on a large optical catalog, drawn 
objectively from a three dimensional map of the large scale 
structure, is essential for understanding the physical prop- 
erties of galaxy groups. Here we construct the first such 
catalog. Our goals are (1) to investigate similarity break- 
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ing in the fundamental scaling laws of systems of galaxies, 
(2) to make the first calculation of the X-ray selection func- 
tion of galaxy groups, and (3) to place firm limits on the 
fraction of optically selected groups that are bound. 

One important example of similarity breaking is the 
relationship between the X-ray luminosity Lx and the 
average plasma temperature T. Ponman et al. (1996) 
show that the Lx — T relation is quite steep for com- 
pact groups, with Lx oc T 5 , whereas for rich clusters 
Lx oc T 3 . This result is consistent with a "preheating" 
scenario where winds from supernovae in galaxies under- 
going the starburst phase leave their mark on the poorest 
systems. Such winds would deplete the intragroup plasma 
(Davis, Mulchaey, & Mushotzky 1999; Hwang et al. 1999), 
raise the gas entropy relative to the gravitational collapse 
value (Ponman, Cannon, & Navarro 1999), and preferen- 
tially dim the systems with the lowest temperatures (Cav- 
aliere, Menci, & Tozzi 1997). 

Systems in hydrostatic equilibrium should have T oc cr 2 , 
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where o~ p is the velocity dispersion of the dark matter halo 
in which the galaxies are embedded. Thus one might ex- 
pect that the Lx — <Jp and the Lx — T relations for groups 
of galaxies steepen in a similar manner. Here we show that 
quite the opposite is true. The groups with the smallest 
velocity dispersions are in fact overluminous compared to 
the Lx oc cr 4 law valid for higher velocity dispersion sys- 
tems. Thus the similarity breaking in the Lx — T law is 
apparently incommensurate with the break in the Lx — &p 
relation. We discuss several plausible explanations for this 
lack of concordance. 

The paper is organized as follows. After constructing 
the catalog (§2), we examine the Lx — o> relation (§3), 
calculate the selection function (§4), discuss the L x — cr p 
flattening (§5), and summarize our findings (§6). We call 
our groups the ROSAT All-Sky Survey— Center for As- 
trophysics Loose Systems, or RASSCALS. 

2. DATA 
2.1. Optical Group Selection 

We extract the optical group catalog for the RASSCALS 
study from two complete redshift surveys. Our catalog in- 
cludes a wide variety of systems, from groups with only 5 
members to the Coma cluster 5 . 

The Center For Astrophysics Redshift Survey (Geller & 
Huchra 1989; Huchra et al. 1990; Huchra, Geller, & Cor- 
win 1995; CfA) and the Southern Sky Redshift Survey 
(Da Costa et al. 1994; Da Costa et al. 1998), both com- 
plete to a limiting Zwicky magnitude m z ~ 15.5, serve as 
sources for the RASSCALS. The portion of the surveys we 
use covers one fourth of the sky in separate sections de- 
scribed in Table |]. We transform the redshifts to the Local 
Group frame (Acz — 300 sin I cos 6), and correct them for 
infall toward the center of the Virgo cluster (300 km s _1 
towards a 20 oo = 12 h 31.2m, <5 2 ooo = 12°2.54'). 

We use the two-parameter friends-of-friends algorithm 
(FOFA) to construct the optical catalog. Huchra & Geller 
(1982) first described the FOFA for use with redshift sur- 
veys, and Ramella, Pisani, & Geller (1997) applied it to 
the NRG data. The FOFA is a three-dimensional algo- 
rithm which identifies regions with a galaxy overdensity 
Sp/p greater than some specified threshold. A second 
fiducial parameter, Vq, rejects galaxies in the overdense 
region which are too far removed in velocity space from 
their nearest neighbor. The N-body simulations of Fred- 
eric (1995) and Diaferio (1999) show that the Huchra & 
Geller (1982) detection method misses few real systems, 
at the cost of including some spurious ones. We apply 
the FOFA to the combined NRG, SRG, and SS2 redshift 
surveys with 5p/p = 80. 

The RASSCALS optical catalog contains 260 systems 
with n > 5 members and 3000 km s _1 < cz < 12000 
km s _1 . The low velocity cutoff rejects systems that cover 
a large area on the sky and thus may be affected by the 
Local Supercluster. The median recession velocity for the 
systems is 7000 km s^ 1 ; the effects of cosmology and evo- 
lution are negligible throughout the sample. Table || lists 
the individual groups and their properties. Figures [| [| 
show the sky positions of the member galaxies for the sys- 

5 In a previous work (Mahdavi et al. 1999), we referred to the Cente 
from the RASSCALS, which was to be the X-ray catalog. We no longi 
as the RASSCALS. 



terns with statistically significant extended X-ray emission 
in the RASS. 

To compare the membership of groups which have dif- 
ferent redshifts we also compute nn, the number of group 
members brighter than an absolute magnitude M z = —17, 
corresponding to m z — 15.5 for a group at cz = 3200 
km s . To calculate nn, we assume that the galaxies in 
groups have the same luminosity function as the Center 
for Astrophysics Redshift Survey (Marzke et al. 1994), re- 
convolved with the magnitude errors. The resulting distri- 
bution is well-represented by a Schechter (1976) function 
with a characteristic absolute magnitude Af* = —19.1 and 
a faint-end slope a = — 1. Table @ lists nn, which has a 
median value of 44. 

2.2. X-Ray Field Selection 

For every system in the RASSCALS optical catalog, we 
obtain X-ray data from a newly processed version of the 
ROSAT All-Sky Survey (Voges et al. 1999), which cor- 
rects effects leading to a low detection rate in the original 
reduction. 

We first assign each system a seven-character name, be- 
ginning with "NRG," "SRG," or "SS2," followed by "b" 
or "s" (specifying the angular size of the system as "big," 
with cz < 8500 km s _1 or "small," with cz > 8500 km s _1 , 
respectively), followed by a three-digit number. 

For each "small" system we extract a square field mea- 
suring 2° x 2° from the RASS; for the "big" systems we 
extract a 3.5° x 3.5° square. The fields are centered at 
the mean RA and DEC of the galaxies; every field is at 
least large enough to include a circle with a projected ra- 
dius of I/i^qq Mpc around the optical center of the system 
it contains. We use photons in the 0.5-2.0 keV hard en- 
ergy band of the Position-Sensitive Proportional Counter 
(PSPC channels 52-201). 

2.3. Detection Algorithm 

The X-ray detection algorithm consists of four steps: 
measurement of the background, source identification, de- 
contamination, and measurement of the source flux. 

We determine the mean background by temporarily re- 
binning the exposure- weighted ROSAT field into an image 
with 15' pixels. We clean the image of all fluctuations with 
an iterative, 2.5er clipping algorithm. The adopted back- 
ground is then the average of the remaining pixel values. 

To estimate the probability that a given group is an 
X-ray source, we use an optical galaxy position template 
(GPT). Mahdavi et al. (1997), who search the RASS for 
X-ray emission from a small subset of our sample, describe 
this method in greater detail. The GPT is defined as the 
union of all projected d = 0.2/i^qq Mpc regions around 
the group members, excluding any galaxies isolated by 
more than d from the rest of the group. We count the 
X-ray photons within the GPT and evaluate the probabil- 
ity that they are drawn from the background distribution. 
All groups that have a detection significance greater than 
2.5<7 progress to the next step. 

We identify the emission peak which coincides most 
closely with the optical center of the group as its X-ray 
counterpart, and calculate the X-ray position of the group 

for Astrophysics— SSRS2 Optical Catalog (CSOC) as a distinct entity 
make that distinction, and refer to the X-ray/optical catalog simply 
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with the intensity-weighted first moment of the pixel val- 
ues. Using standard maximum likelihood techniques, we 
identify contaminating X-ray point sources over the entire 
field. We remove these sources by excising a ring of radius 
3' , roughly three times the full width at half maximum of 
the ROSAT PSPC point spread function (PSF). Unrelated 
extended sources often contaminate the group emission; we 
use a suitably larger aperture to remove them. We have 
examined publicly available ROSAT High-Resolution Im- 
ager (HRI) observations of a few groups, and find that our 
RASS decontamination procedure is satisfactory. 

To reject groups with entirely pointlike X-ray emis- 
sion, we calculate N(R), the cumulative distribution of 
the ROSAT counts. We use use the Kolmogorov-Smirnov 
(KS) Test to compare the shape of the emission peak with 
that of the PSPC PSF combined with the background. 
We take sources with Pks < 0-05 as inconsistent with the 
PSF. 

Finally, we convert the PSPC count rate into Lx(R), 
the 0.1-2.4 keV X-ray luminosity contained within a ring 
of projected radius R. The Appendix describes the flux 
conversion procedure in detail. 




Fig. 1.— Properties of the NOCORE estimator: (a) for 
/3-models with /3 = 0.65 (solid line) and f) = 1 (dashed line); (b) R^, 
the radius where £,{R) is minimum, as a function of /3. 



2.4. Core Radius Estimation 

Here we describe a procedure for identifying a physical 
scale for the X-ray emission. There is a great deal of ev- 
idence that the emissivity profiles of clusters of galaxies 
exhibit a characteristic scale, or core radius, r c (e.g. Jones 
& Forman 1984; Mohr, Mathiesen, & Evrard 1999). For 
example, the "/3-model" emissivity profile frequently used 
to fit observations of dynamically relaxed systems, 



e(r) oc II 



-3/3 



(1) 



files (tar 1 for r <C r c ) in clusters with cooling flows 
(Thomas 1998). 

The usual method for measuring r c from X-ray observa- 
tions consists of projecting e(r) along one dimension, and 
fitting the resulting surface brightness profile to the data. 
This approach has the disadvantage that the resulting es- 
timate of r c is model-dependent, and is often strongly cor- 
related with the slope parameter /3, even with very high 
quality data (e.g. Jones & Forman 1984; Neumann & Ar- 
naud 1999). Furthermore, its application to the RASS- 
CALS is limited, because the small number of counts and 
the relatively large uncertainty in the background make it 
difficult to reconstruct accurate surface brightness profiles 
for all but the brightest systems. 

We therefore use the Nonparametric Core Radius Esti- 
mator (NOCORE; Mahdavi 2000) to avoid the core fitting 
procedure and its associated uncertainties. NOCORE is 
model-independent; it does not require an estimation of 
the background, and it relies on the properties of the inte- 
grated emission profile, rather than the differential profile, 
to estimate the core radius. Its only assumption is the 
constancy of the background level at the position of the 
object of interest. 

Now we outline the procedure. Consider the measured 
count rate within an annulus R from the X-ray center of 
the group: it consists of the emission of the group itself, 
S(R), plus the constant background count rate per unit 
area, B: 



N{R) = S{R) 



itR 2 B. 



(2) 



The fundamental basis of NOCORE is the observation 
that the quantity N(R) — k 2 N(R/k), where k is a number 
greater than 1, is completely independent of the constant 
background. Formally, we define the NOCORE radius as 
the radius where the function 



_ N(R) - AN{R/2) 
= R 



(3) 



is nearly constant for physical radii r <C r c , and scales as 
r -6/3 f or r Tc There is also evidence for cuspy pro- 



has a global minimum. The division by R is necessary to 
obtain a detectable minimum. 

Figure [j]a shows £(i?) for theoretical j3- models with 
(3 = 0.65 and f3 — 1. The minimum, is well-defined 
in both cases. The location of R^ as a function of f3 ap- 
pears in Figure [lb. We have carried out numerical tests 
of the method, adding Poisson noise and a background to 
a variety of /3-models, to verify that we recover the ap- 
propriate core radius without bias. When applying the 
method to the observations, we use bootstrap resampling 
to determine 68% confidence intervals on R^. 

As long as there is a characteristic scale in the emissiv- 
ity profile of a system, NOCORE will find it. The func- 
tion £(i?) has a well-defined minimum even in cases when 
the /3-model is not a good description of the emissivity, 
for example systems with a cooling flow. If there is more 
than one characteristic scale in the profile — if the emissiv- 
ity has features at several radii, because of substructure 
in the cluster, for example — then £(i?) has more than one 
minimum. If the profile is a pure power law, NOCORE 
shows no core; £(i?) is then a monotonically decreasing or 
increasing function of R. 

We use the radius at which £(i?) is minimum as a mea- 
sure of the physical scale of the X-ray emission. 
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Fig. 2. — The RASSCALS. Sky maps are 1.5/i 10 q Mpc wide, and include the name of each group, its mean recession velocity, velocity 
dispersion, and log X-Ray luminosity within 0.5/iJ~qq Mpc. The solid X-ray emission contours begin at ltr above the background and increase 
by a factor of 2. Dotted regions represent sources excised by the decontamination algorithm. Circles represent early-type galaxies, and the 
elongated symbols represent late-type galaxies. The size of the symbols is proportional to their apparent Zwicky magnitude; the largest 
symbol corresponds to mz = 12, and the smallest corresponds to mz = 15.5. Also shown are the background-subtracted cumulative X-ray 
luminosity profile, and the NOCORE estimator £;(/?). The 68% confidence interval of the NOCORE radius appears as a horizontal bar. 
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Fig. 3. — See caption for Figure^. Figures for the remainder of the X-ray detected RASSCALS will appear in the published version of this 
article, in The Astrophysical Journal. 



6 



Mahdavi et al. 



2.5. General Properties of the Final Catalog 

Table || lists the projected velocity dispersion, a p , the 
0.1-2.4 keV X-ray luminosity, and the NOCORE radius 
i?f of the detected groups. There are 61 detections, of 
which two (NRGs372 and NRGs392) are bright X-ray clus- 
ters (Abell 2147 and 2199) that the friends-of-friends algo- 
rithm has broken up into pieces. We count these clusters as 
detections but we do not calculate luminosities for them. 
Figures show galaxy positions and X-ray emission con- 
tours for the detected groups. 

Because the diffuse X-ray emission is generally a marker 
of gas held in a gravitational potential, these systems 
are probably bound configurations. There is, however, a 
chance that the X-ray emission might be due to projection 
along an unbound filament in the large scale structure of 
the universe (Hernquist, Katz, & Weinberg 1995). The 
galaxies projected along this line of sight might also have 
similar redshifts without being bound. However, deeper 
redshift surveys in the fields of X-ray emitting RASSCALS 
show that a number of these systems have velocity disper- 
sion profiles o p (R) that decline as a function of projected 
distance from the group center; these profiles are consis- 
tent with the expectation for a relaxed dynamical system 
(Mahdavi et al. 1999). 

In Figures |]^3] we also show the cumulative luminosity 
profile and the NOCORE estimator £(i?). In several cases 
£(i?) appears to have several local minima in addition to 
the global minimum. The local minima are almost always 
due to Poisson noise and deviations from spherical symme- 
try in the structure of the gas. The confidence intervals on 
take these fluctuations into account: when the fluctua- 
tions in £(i?) dominate its shape, the error in is large. 
But when £(i?) has a well-defined global minimum and the 
fluctuations are small, i?j is relatively well determined. 

3. L x - op RELATION 

Here we examine the relationship between the X-ray lu- 
minosity and the projected velocity dispersion. First we 
comment on the link between the Lx — o~ p scaling law, 
which relates X-ray and optical data, and the Lx — T scal- 
ing law, which is internal to X-ray data. Then we describe 
the actual relation. 

3.1. Background on the Scaling Laws 

If the member galaxies trace the total mass distribution 
in a cluster, a simple theoretical calculation (Quintana & 
Melnick 1982) predicts that a spherically symmetric ball of 
gas should have Lx oc pa^T 1 / 2 , where / is the ratio of the 
gas mass to the total mass. A further, common assumption 
is that T oc cr 2 , i.e., that the emission- weighted gas tem- 
perature is proportional to the depth of the gravitational 
potential. These assumptions yield Lx oc f 2 T 2 oc / 2 cr 4 . 

The observed Lx — er p relation for rich clusters is in 
good agreement with the theoretical prediction; Quintana 
& Melnick (1982) and Mulchaey & Zabludoff (1998), for 
example, find slopes consistent with Lx oc a 4 . The empir- 
ical Lx — T, relation, however, is somewhat steeper than 
expected, with most finding Lx oc T 2 - 75 , even after re- 
moving the central cooling flow region (e.g., Markevitch 
1998 and references therein). 

If the discrepancy between the simple theoretical predic- 
tion, Lx oc T 2 , and the observations is real, several effects 



might explain it. It could be that / increases slightly with 
T (David, Jones, & Forman 1995), or that T oc a 1 - 5 , con- 
sistent with a nonisothcrmal, polytropic gas distribution 
(Wu, Fan, & Xu 1998). Finally, preheating of gas in the 
kT ^ 4 keV systems may preferentially dim them, leading 
to a steeper relation. Ponman, Cannon, & Navarro (1999) 
show that this latter possibility is particularly attractive 
because it also accounts for differences in the shapes of 
X-ray surface brightness profiles among kT ^ 4 keV and 
kT ^ 4 keV clusters. Cavaliere et al. (1997) work out the 
Lx — T relation for this scenario, and find that it steepens 
gradually as T declines, with Lx oc T 5 for poor groups , 
L x oc T 3 for 2 keV < kT <, 7 keV systems, and L x oc T 2 
for the hottest clusters. This Lx — T relation fits temper- 
atures and luminosities for a range of systems from poor 
groups to clusters. 

Now, if a single power law describes the scaling of the 
velocity dispersion o~ p with the temperature T, and the 
Cavaliere et al. (1997) preheating model is correct, one 
should observe a similarly steep Lx — cr p relation for poor 
groups. Three different works have attempted a measure- 
ment of the faint end of this relation, with three different 
results. 

1. Ponman et al. (1996) analyze a mixture of pointed 
and RASS observations of a sample of Hickson (1982) 
Compact Groups (HCGs hereafter). They obtain 
Lx oc crp 9±21 for the groups with pointed obser- 
vations. While this result favors a steeper slope, the 
68% confidence interval is quite large: the HCGs 
contain as few as three member galaxies, and hence 
it is very difficult to estimate the correct velocity dis- 
persion. Also, HCGs are often embedded in much 
richer systems (Ramella et al. 1994), and this em- 
bedding may further bias the value of the velocity 
dispersions. 

2. Mulchaey & Zabludoff (1998; MZ98 hereafter) carry 
out deep optical spectroscopy for a more limited 
sample of poor groups with pointed ROSAT observa- 
tions. Because they obtain « 30 members per group, 
their derived velocity dispersions should be more re- 
liable than those of Ponman et al. (1996). They ob- 
tain Lx oc (7 4 - 3±0 - 4 for a combined sample of groups 
and clusters. 

3. Mahdavi et al. (1997) use our method to examine 
ROSAT data for a small but statistically complete 
subset of the RASSCALS. They do not, however, ex- 
cise emission from individual galaxies; furthermore, 
they assume a constant plasma temperature kT = 1 
keV, rather than leaving T free to vary as we do here. 
They obtain Lx oc <7p' 56±a25 , much shallower than 
either Ponman et al. (1996) or MZ98. 

In summary, Ponman et al. (1996) find an Lx — <J P rela- 
tion consistent with the simplest predictions of preheating 
models; MZ98 derive a relation that is consistent with the 
standard picture with no preheating; and Mahdavi et al. 
(1997) find that the faint-end slope is much shallower than 
the prediction of either of the two scenarios. 

We now consider the Lx — o~ p relation for the com- 
plete set of RASSCALS. Our procedure differs from that of 
Mahdavi et al. (1997), because we remove contaminating 
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sources whenever they are detectable, model the plasma 
temperature, and use an updated version of the RASS. 

To compare the RASSCALS Lx — cr p relation with that 
of richer systems, we take cluster X-ray luminosities from 
the paper by Markevitch (1998), where cooling flows are 
removed from the analysis. We use only clusters which 
have velocity dispersion listed in Fadda et al. (1996), who 
consider systems with > 30 measured redshifts. Table || 
lists these data. Figure || shows the combined cluster- 
RASSCALS data. The Lx — cr p seems to flatten as the 
luminosity decreases. 

3.2. Details of the Fitting Procedure 

To place a quantitative constraint on the degree of flat- 
tening, we fit a broken power law of the form 



log — = s {o- pi o- k ) log— ; 

. . / si if Op < cT fe 
s(er„, a k ) = < . t 
y P ' K > 1 32 if cr p > a k 



(4) 
(5) 



Here si and S2 are the faint-end and bright-end slopes, 
respectively, and (a k , L k ) is the position of the knee of the 
power law. We then minimize a merit function appropriate 
for data with error in two coordinates (Press et al. 1995, 
§15.3), 



[\og(Lj/L k ) - s{ai,a k )\og (ajj 'o- k )] 2 
—[ (A log L^ 2 + s((T j; ,(T fc ) 2 (A logo-;) 2 



(6) 



where {o~i,Li) are the measurements, with errors 
(A<7 i; ALj). To minimize the % 2 we apply the following 
procedure. 

1. First, we fit a single power law by forcing s\ = S2 
and log Cfc = 0, and applying the Press et al. (1995, 
§15.3) package. The result appears as the dashed line 
in Figure ||. We call this best-fit power law slope Sq. 

2. Next, we vary the position of the knee of the power 
law over a 50 x 50 grid with bounds \oga k = [2,3] 
and logifc = [42,44]. At each point in the grid, we 
minimize the x 2 over si and S2, using the Fletcher- 
Reeves-Polak-Ribiere algorithm, which makes use of 
gradient information (Press et al. 1995, §10.6). We 
start the minimization algorithm with s\ = §2 = sq, 
and require < si < 10 and < S2 < 10 as priors. 
This procedure yields a function Xmin( a k, L k ). 

3. Finally, we minimize Xmin(. a k, L k ) to obtain the po- 
sition of the knee of the power law and the best-fit 
slopes associated with it. The results of the fit ap- 
pear in Figure ^. 

We also try more robust estimators, such as the BCES 
bisector (Akritas & Bershady 1996). In general, these es- 
timators are in good agreement with the results of the 
X 2 fits; however, Akritas & Bershady (1996) do not pro- 
vide a mechanism for assessing the quality of the fit, and 
their package does not allow for the calculation of joint 
two-dimensional confidence intervals. We therefore focus 
on the x 2 statistic. Below we also consider how the fit 
changes with the inclusion of the 199 upper limits. 



3.3. A Broken Power Law is the Best Fit 

Our data unambiguously favor a broken power law over 
a single power law for the Lx — o v relation. The confi- 
dence contours in Figure || show that the faint-end slope 
and the bright-end slope are different at better than the 
99.7% confidence interval. Furthermore, the scatter in the 
Lx — o~p relation is actually reduced by fitting a broken 
power law instead of a single power law. 

The faint-end slope, s\ — 0.37 ± 0.3, is even shallower 
than the earlier finding of Mahdavi et al. (1997), si = 
1.56 ± 0.25, for their fit to 9 low-luminosity RASSCALS. 
The shallowness of our faint-end slope is remarkable be- 
cause, unlike Mahdavi et al. (1997), we remove sources 
of individual emission whenever possible, and model the 
plasma temperature without fixing it at a particular value. 
The bright-end slope, S2 = 4.02 ±0.1, on the other hand, 
is consistent with MZ98, whose Lx — o v depends mainly 
on rich clusters; their single power law fit has a slope 
4.29±0.37. A slight discrepancy is to be expected, because 
MZ98 use bolometric X-ray luminosities, and we measure 
the luminosities in the 0.1-2.4 keV spectral range. How- 
ever, this discrepancy should cause only an As2 ~ 0.4 
offset in the slope for clusters with Lx 10 43 ergs s . 
Systems with Lx Ss 10 43 ergs s _1 should have bolometric 
luminosities comparable to their 0.1-2.4 keV luminosities. 

We stress that our fitting procedure in no way favors the 
shallower slope: we begin the x 2 minimization by setting 
both slopes equal to the best-fit single power law. Also, a 
broken power law is the best fit even if we exclude the low- 
est velocity dispersion group, SRGb075, from the fit. Do- 
ing so, we would obtain si = 1.39±0.5 and S2 — 3.99±0.3. 

Finally, we consider whether including the upper 
limits in the fit changes the derived slopes. For 
this task we obtain the Astronomy Survival Analysis 
Pack age (ASURV; Lavalley, Isobe, fc F eigelson 1992) 
from http://www.astro.psu.edu/statcodes. ASURV im- 
plements the methods described in Isobe, Feigelson, & 
Nelson (1986) for regression of data which includes both 
detections and upper limits. Furthermore, ASURV allows 
for an intrinsic scatter in the relation. 

For the 155 objects with a p < 340 km s _1 , we find that 
the best-fit slope is 1.38 ± 0.4; for the 128 objects with 
cr p > 340 km s -1 , it is 5.37 ± 0.5. Thus the inclusion 
of upper limits does not bring the two slopes closer to 
each other; if anything, it strengthens our claim that the 
Lx — o~p relation is best described by a broken power law. 



4. DETECTION STATISTICS 

We now examine the statistical properties of the catalog 
in greater detail. We seek a deeper understanding of the 
X-ray selection function of the RASSCALS. A useful tool 
for this purpose is the number distribution of a set of mea- 
surements x, which we label N(x). For example, we call 
the number distribution of the group velocity dispersion 
N{<T p )] the number distribution of the distance-corrected 
group membership is N(ni7). 

Although the traditional estimator of the number dis- 
tribution is the histogram, we compute N{x) using the 
DEDICA algorithm (Pisani 1993). DEDICA makes use of 
Gaussian kernels to arrive at a maximum-likelihood esti- 
mate of the number distribution. The resulting smooth 
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Fig. 4. — The Lx — Cp relation. The dashed line shows the single power law fit, and the solid line shows the broken power law. The 
bottom panel shows the joint 68.3%,and 99.7% confidence intervals for the slopes of the broken power law. The undetected RASSCALS are 
not considered in the fit, but in §3.3 we consider the effect of including the upper limits. 
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function, N(x) is more useful than a histogram, because 
N(x) is nonparametric, and any structure within N(x) is 
statistically significant. We normalize N(x) so that the 
number of groups with x\ < x < x^ is given by, 



N(x)dx. 



(7) 



4.1. An Abundance of Low o~ p Systems 

Figure ||a shows N(a p ) and N(nn) separately for all 
RASSCALS and for those with significant X-ray emission. 
Interestingly, N(a p ) for all 260 groups is double peaked; 
there are 102 RASSCALS (39%) with a p < 250 km s" 1 . 

The abundance of these low a p systems is puzzling con- 
sidering that many of them probably contain unrelated 
galaxies ("interlopers") (Frederic 1995), and that these 
interlopers typically lead to an overestimate, not an un- 
derestimate, of the velocity dispersion. Several plausible 
explanations for their frequent occurrence exist. 

1. It may be that the groups with low a p are unbound, 
chance projections along the line of sight. However, 
this situation is highly unlikely for systems drawn 
from a complete redshift survey. Chance superposi- 
tions in such a survey in fact have a larger mean ve- 
locity dispersion than the bound groups do (Ramella 
et al. 1997). 

2. The groups might be pieces of sheets or filaments in 
the large scale distribution of matter. A collapsing 
sheet of galaxies which is removed from the Hubble 
flow, and which is perpendicular to the line of sight, 
might look like a group to the friends-of-friends al- 
gorithm. 

3. The velocity dispersion a p might not be related to 
the mass distribution in a straightforward manner. 
For example, Mahdavi et al. (1999) find that galaxy 
orbits in a subsample of the RASSCALS have a sig- 
nificant mean radial anisotropy; and Diaferio (2000) 
uses N-body simulations to show that systems with 
cr p ^ 300 km s _1 have galaxy velocity dispersions 
that are uncorrelated with the total group mass. 



4.2. Detection Ef 



iciency 



Here we use the groups we have detected as a basis for 
estimating (1) the number of RASSCALS with X-ray emis- 
sion too faint to be observable by ROSAT, and (2) the 
number of RASSCALS with no X-ray emission, some of 
which might be unbound superpositions. 

We begin by assuming that all the RASSCALS emit X- 
rays according to an empirical relationship between Lx 
and a p . We compute the number of groups we expect to 
detect as a function of a p , and compare the theoretical 
detection probability with the true detection efficiency. 

Suppose that all the RASSCALS emit X-rays according 
to a power law relationship between Lx and cr p , 



log L x = s log a p + b. 



(8) 



If the local flux detection threshold for a group at redshift 
z is Fq, it will be detectable if a p > erg, where 

log[47rF (l + z) 2 c 2 z 2 ] -b 
log cr = - - (9) 



The theoretical probability that the group will be detected 
is then 

/>oo 

Pth(cr p ) = / p(l0gCT p )dl0gCTp, (10) 

J log (70 

where p(\og a p ) is the probability distribution function of 
log cr p . We calculate P th (a p ) for all 260 RASSCALS, tak- 
ing s — 4.02 and b = 32.19 from the single power law 
determined in §||. We approximate log a p as a Gaussian 
with a standard deviation equal to 1.3 times the uncer- 
tainty given in Table || for each group. Multiplying the 
error in a p by 1.3 is a way of spreading the uncertainty in 
the Lx — o- p relation directly into Pth(op). The resulting 
average theoretical probability of detecting a group with 
velocity dispersion a p is well approximated by 



Pth(o- p ) 



1 1 



rf 



(log 



250 km s" 



(11) 



where erf (a;) is the error function. 

Figure |jb shows the observed and the theoretical detec- 
tion probabilities. The solid line represents the fraction 
of the RASSCALS we actually detect, P o bs(0>), and the 
short dashed line shows Pth(o>), the fraction of the RASS- 
CALS we should detect given Lx oc er 4 . The quotient, 



, , s _ Pobs(o-p) 



(12) 



appears as the long dashed line. fx{o~ p ) represents the 
fraction of groups that should have extended X-ray emis- 
sion in order that we detect our set of 61 RASSCALS. 
Remarkably, fx is a nearly constant 40% for a p > 150 
km s _1 , and rises steeply for a p < 150 km s _1 . The scat- 
ter around the theoretical probability Pth(tfp) introduces 
a 30% uncertainty in the breaking point a p = 150 km s _1 , 
but does not affect the result that fx ~ 40% ± 8% above 
the breaking point. 

Thus Figure |^b shows that we detect many fewer sys- 
tems overall than expected from the raw Lx oc a p relation. 
To match our observed detection efficiency, only 40% of 
groups with a p > 150 km s -1 must have extended X-ray 
emission. The probability fx that a group contains X- 
ray emitting gas does not seem to increase with the group 
velocity dispersion. 

On the other hand, the detection of the a p < 150 km s _1 
groups exceeds the expectation from Lx oc cr 4 . The theo- 
retical probability of detecting any of these low-cr p groups 
is near 0, and yet we detect 6% of them. A flattening of the 
true Lx — cr p relation for low velocity dispersion systems, 
of the kind we discuss in §0, resolves the discrepancy. 

The result that only 40% of the a p > 150 km s" 1 RASS- 
CALS should emit X-rays has an interesting interpretation 
when combined with the predictions of N-body (Frederic 
1995, Diaferio 1999) and geometric (Ramella et al. 1997) 
simulations of the local large-scale structure. These sim- 
ulations suggest that <^ 80% of groups with n > 5 mem- 
bers drawn from a complete redshift survey should be real, 
bound systems. If indeed > 80% of the RASSCALS are 
bound, and our simulations are correct, then at least half 
the bound groups must possess a negligible amount of ex- 
tended X-ray emission. 

The X-ray data impose a lower limit of 40%, and the 
simulations impose an upper limit of 80%, on the fraction 
of RASSCALS that are real, bound systems of galaxies. 
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Fig. 5. — Number distributions of a v and ni?: (a) distribution of a v for all groups (dashed line) and for the Q = A, B sample (solid line); 
(b) the observed detection efficiency as a function of a p (solid line), the expected detection probability assuming Lj oc a 4 (short dashed 
line), and their quotient (long dashed line). Similar analysis for ni7, the number of system members brighter than an absolute magnitude 
M z = —17, appears in (c) and (d). The dotted horizontal line shows the mean detection efficiency for the RASSCALS (23%). 
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5. DISCUSSION 

The flattening of the Lx — o p relation for systems with 
o p ^ 340 km s _1 is now well established, not just by our 
study, but by pointed iJOS'^lT'observations of an indepen- 
dent sample of 24 groups (Helson & Ponman 2000) . This 
similarity breaking is particularly striking, because it is in 
conflict with the Lx ~ T relation for systems of galaxies, 
which actually steepens as the temperature drops (Metzler 
& Evrard 1994; Ponman et al. 1996). 

It is difficult to dismiss the Lx — cr p flattening by claim- 
ing that the velocity dispersion of the discrepant groups is 
biased towards lower values. Most groups drawn from red- 
shift surveys contain unrelated galaxies which tend to bias 
<r p towards larger values (Frederic 1995; Diaferio 1999). 
It is also no longer possible to argue (e.g. Mulchaey & 
Zabludoff 1998) that the flattening is due to a failure to 
remove detectable contamination. We, as well as Helson & 
Ponman (2000), remove such contamination to the extent 
allowed by the data. 

The shallow Lx — <J P slope may be explainable through 
the "mixed-emission" scenario proposed by Dell' Antonio 
et al. (1994). It is possible that a number of galaxies with 
faint, X-ray emitting ISMs are embedded within the in- 
tragroup medium. These individually emitting galaxies 
could contribute significantly to the total luminosity and 
place it above the virial value. A large fraction of such 
emission would be neither directly detectable nor remov- 
able, appearing instead as fluctuations in excess of Poisson 
and instrumental noise superposed on the central emission 
peak (Soltan & Fabricant 1990). Further verification of 
the mixed emission scenario thus depends on higher qual- 
ity observations of the lowest velocity dispersion systems 
with the Chandra or XMM missions. 

However, we can investigate whether the break in the 
Lx — cr p relation is linked to other physical properties of 
the RASSCALS. One possibility is that the excess emis- 
sion is characteristic of the dynamically youngest groups, 
those perhaps still in the process of formation. An in- 
dicator of such a dynamical state might be the fraction 
of spiral member galaxies, / sp . If the dominant process 
for the formation of elliptical galaxies in groups is galaxy- 
galaxy interaction, then one might expect a system with 
/ sp ss to be much more evolved than a group composed 
mainly of spiral galaxies. 

Figure ||a shows a weak correlation between the X-ray 
luminosity and the spiral fraction. The correlation be- 
tween the two quantities is barely significant (Kendall's 
r = —0.109, with P — 0.22, a 1-cr result). It is notewor- 
thy that no system with / sp ; > 0.5 is more luminous than 

5^ioo x 10 42 er S s s_1 — groups that are spiral-dominated 
tend to have below average X-ray luminosities. 

However, closer inspection reveals that the spiral frac- 
tion is not related to the Lx — <J P flattening. SRGb075, 
the X-ray emitting group with the smallest velocity dis- 
persion (tip = 60 km s _1 ) has / sp i = 0.2. The group with 
the next smallest a p , SS2b293, has / sp = 0.33, and the fol- 
lowing group, NRGb045, has / sp = 0.2. Although a trend 
relating / sp and Lx probably exists, the RASSCALS that 
are responsible for the flattening of the Lx — <? P relation 
have spiral fractions comparable to those of higher velocity 
dispersion groups. 

Another possible indicator of the dynamical age of a sys- 



tem of galaxies is its crossing time. Groups where galaxies 
have completed many orbits might be closer to dynamical 
equilibrium than those where the galaxies have made only 
a few crossings. We note, however, that if the accretion of 
external galaxies plays a significant role in the evolution of 
a group, it may not reach dynamical equilibrium even after 
many crossing times have passed (Diaferio et al. 1993). 

The crossing time of a system in units of the Hubble time 
is roughly t c — RHq/(t Pi where R is the characteristic size. 
Because many of the groups in our sample have fewer than 
9 members, computing R from the optical data is likely to 
lead to large errors in t c . Instead, we use the NOCORE 
radius (§2.4) to estimate the crossing time, t c = R^Hq/<j p . 
The NOCORE radius provides a characteristic scale for 
the X-ray emission, and hence for the gravitational poten- 
tial of each group. 

There is a significant correlation between Lx and t c in 
Figure |b (Kendall's t = -0.228, P = 0.01, a nearly 3-cr 
result). Of course, this effect follows directly from the re- 
lationship between Lx and a p (which exhibit a 10- ct corre- 
lation). Because R^ is uncorrelated with er p , the inclusion 
of i?£ increases the scatter. 

However, the Lx — t c comparison does reveal an in- 
teresting property of the groups which contribute to the 
flattening of the Lx — <J p correlation. These groups have 
t c > 0.3HQ 1 ; they have longer crossing times than the 
groups in the steeper, Lx oc cr 4 portion of the relation. 
Thus we have an indication that the X-ray overluminous 
groups are also the ones where the crossing time is a large 
fraction of the Hubble time. An explanation of this result 
in terms of the dynamical histories of the 1ow-ct p groups 
awaits a much deeper optical and X-ray probe of their 
structure. 



6. CONCLUSION 

The RASSCALS are the largest extant combined X-ray 
and optical catalog of galaxy groups. We draw the systems 
from two redshift surveys that have a limiting magnitude 
of m z = 15.5 and cover tt ster of the sky. There are 260 
systems, of which 23% have statistically significant X-ray 
emission in the ROSAT All-Sky Survey after we remove 
contamination from unrelated sources. We include a cat- 
alog of the systems. 

We calculate the X-ray selection function for our sam- 
ple. The behavior of the function implies that only 40% 
of the RASSCALS are intrinsically X-ray luminous. The 
remaining w 60% of the RASSCALS are either chance su- 
perpositions, or bound systems devoid of hot gas. 

We examine the relationship between the X-ray lumi- 
nosity Lx and the velocity dispersion a p for the 59 high- 
quality RASSCALS and a representative sample 25 of rich 
clusters not internal to our data. The best fit relation 
is a broken power law with Lx oc ,j0- 37±0 - 3 for a p < 340 

km s _1 , and L x oc o^ 9 * - 1 for a p > 340 km s" 1 . Whether 
we include the upper limits in our analysis, or assume a 
dominant intrinsic scatter in the relation, a broken power 
law with a shallow faint-end slope is still a better fit than 
a single power law. 

Stressing that we have been careful to remove contami- 
nation from individual galaxies and unrelated sources, we 
conclude that the flattening in the Lx — o p relation for 
groups of galaxies is a physical effect. A potential mecha- 
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nism for the excess luminosity of the faintest systems is the 
"mixed emission" scenario (Dell' Antonio et al. 1994): the 
emission from the intragroup plasma may be irrecoverably 
contaminated by a superposition of diffuse X-ray sources 
corresponding to the hot interstellar medium of the mem- 
ber galaxies. A final explanation of the flattening of the 
Lx — a relation must focus on the detailed X-ray and opti- 
cal structure of the groups with small velocity dispersions 
(<7 P < 150 km s" 1 ). 

We plan to calculate the X-ray luminosity function of 



the RASSCALS soon. Deep optical spectroscopy of these 
systems is already underway, and the first results appear 
in Mahdavi et al. (1999). 
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Italian Space Agency (M. R.). 



APPENDIX 



LUMINOSITY, FLUX AND TEMPERATURE CALIBRATION 

Here we describe the procedure we use to convert the ROSAT PSPC count rate into a 0.1-2.4 keV X-ray luminosity. Our 
procedure does not require fixing or guessing the plasma temperature. Instead, we fold the uncertainty in the temperature 
directly into the derived luminosities. 

The decontaminated source count rate within a ring of projected radius R is 



S = 



ttR 2 



TlR? ^Llcan 



E 



Ei 



-ttR 2 B, 



(Al) 



where A c \ can is the area of the portion of the ring removed during the decontamination process, Ni is the total count rate 
within pixel i, Ei is the exposure time within pixel i, B is the average background count rate per unit area on the sky, and 
the summation is over all pixels within the ring. Pixels which fall partially inside the ring are appropriately subdivided. 
The error in the source count rate, as, is given by 



4 




irR 2 a B ) 



(A2) 



where <jb is the uncertainty in the background. 

Because all our systems have redshift z < 0.04, the 0.1-2.4 keV X-ray luminosity, Lx, is 



(A3) 



The 0.1-2.4 keV flux, F from the GPT is then 

F = C(N H ,T)S. (A4) 

Here C(Nh,T) is a function suited to the ROSAT PSPC Survey Mode instrumental setup which converts the 0.5-2.0 
keV count rate to the appropriate 0.1-2.4 keV flux from a Raymond & Smith (1977) spectrum with the abundance fixed 
at 30% of the solar value. C(Nh, T) depends on Nh, the total hydrogen column density along the line of sight, which we 
compute using the results of Dickey & Lockman (1990), and the emission-weighted plasma temperature, T. 

We cannot accurately determine T independently of F from the RASS data. However, once Nh is fixed, C{Nh,T) 
varies only 15%-20% for 0.3 keV < kT < 10 keV. We therefore fold this uncertainty in T into our calculation of the flux. 

If pc{C) is the probability distribution function (PDF) of C(Nh , T), and ps(S) is the PDF of the source count rate S, 
then the PDF of the flux is (Lupton 1993, pp. 9-10) 



Pf(F) = / 
Jo 



Pc(C)p s (F/C) 



dC_ 

~~c 



(A5) 



If the PDF of the group's emission- weighted temperature is pr(T), then, by the law of transformation of probabilities, 



Pc{C)=p T {T) 



dT 
dC 



Approximating ps(S) as a Gaussian distribution with mean S and standard deviation as(S), we obtain 

2" 



Pf(F) 



1 



Pt(T) 
2na s Jo C{N H ,T) 



exp 



1 ( F/C(N H ,T)-S 

2 ^ a s 



dT. 



(A6) 



(A7) 



We take pt{T) to be (9.7 keV) -1 over the range 0.3-10 keV, and zero everywhere else. We have also tried a more 
sophisticated approach, with p(T) proportional to the observed temperature function of systems of galaxies (Markevitch 
1998). The difference between the resulting PDF and the constant pt(T) PDF is negligible compared with the error 
introduced by the uncertainty in the temperature function itself. 
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Table 1 

Sky Coverage of the Optical Group Catalog 



Field Name a 2 ooo £2000 



NRG 


8.5 hr 


17 hr 


8.5° 


- 43.5° 


SRG 


21.5 hr - 


3 hr 


-2° 


32° 


SS2 a 


21 hr - 


5 hr 


-40° 


1.5° 


SS2 b 


10 hr - 


15 hr 


-26° 


0° 



a First portion, excluding sections with galactic 
latitude b < 40°. 

a Second portion, excluding sections with galac- 
tic latitude b < 40°. 



The RASSCALS 



Table 2 

The RASSCALS: Basic Properties 



RASSCALS 
ID 


J2000 


<5 a 
J2000 


n 


7117 


cz 
km s _1 


loga p 
km s _1 


logL x b 

Koo er s s_1 


Other 
Identification 


SS2b003 


00:08:55.8 


-37:28:15 


5 


35 


8357+60 


2.14+0.20 


< 42.1 




SRGb061 


00:11:44.7 


+28:21:35 


10 


57 


7855+163 


2.71+0.06 


< 41.9 


PPS058 


SS2b004 


00:14:47.2 


-07:14:16 


5 


11 


5290+65 


2.19+0.14 


< 41.7 




SS2b005 


00:15:31.5 


-24:07:38 


5 


24 


7390+35 


1.86+0.15 


< 42.0 




SRGb062 


00:18:25.2 


+30:04:13 


13 


49 


6811+122 


2.64+0.10 


42.66+0.07 


MGBR 


SRGb063 


00:21:38.4 


+22:24:20 


10 


25 


5665+87 


2.46+0.11 


42.56+0.09 


PPS062 



a For the 59 groups with a listed X-ray luminosity, we report the X-ray centroid; for the others we report the 
mean RA and DEC of the galaxies in the group. 

b Luminosity in the 0.1-2.4 keV band within an aperture of 0. 5/1^00 Mpc. 

c By no means complete. A: Abell cluster; AWM: Albert et al. (1997) Groups; HCG: Hickson Compact Groups 
(Hickson 1982); MGBR: Studied in greater detail in Mahdavi et al. (1999); PPS: Loose Groups in the Perseus- 
Pisces Survey (Trasarti-Battistoni 1998); ZM: Zabludoff & Mulchaey (1998) Poor Groups 



Note. — The complete table will be available in the electronic version of The Astrophysical Journal. The first 
few lines are shown to elucidate form and content. 



Table 3 

The RASSCALS: Detailed Properties 



RASSCALS 


logo-p 


log L x (0.25) 


log L x (0.5) 


logLjr(fle) 


Re 


logi c 


/sp 


ID 


km s _1 




Koo cr s s_1 




h ioo M P C 






SRGb062 


2.64+0.10 


42.41+0.08 


42.66+0.07 


42.27+0.08 


0.15+0.07 


-1.45+0.21 


0.15 


SRGb063 


2.46+0.11 


42.32+0.09 


42.56+0.09 


42.71+0.10 


0.84+0.30 


-0.54+0.19 


0.40 



Note. — Data for the detected groups, excluding the optically spoiled clusters NRGs372 and NRGs392. 
The 0.1-2.4 keV X-ray luminosities are computed within 0. 25/2.700 Mp c and 0.5/i^qq Mpc, as well as within 
the NOCORE radius, R^. The crossing time t c = R^Ho/a p , and / sp is the fraction of group members that 
are spiral galaxies. 

Note. — The complete table will be available in the electronic version of The Astrophysical Journal. The 
first few lines are shown to elucidate form and content. 
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Table 4 
Clusters Data 



l\Tn ttip 


2 


lno" T lV a 

"ioo er § s 


1 no' rj ^ 

km s ^ 


A85 


0.052 


44.28 


2.99±0.02 


A119 


0.044 


43.94 


2.83±0.06 


A399 


0.072 


44.23 


3.05±0.04 


A401 


0.074 


44.46 


3.06±0.03 


A478 


0.088 


44.51 


2.96±0.16 


A754 


0.054 


44.35 


2.82±0.04 


A1651 


0.085 


44.26 


3.00±0.07 


A1736 


0.046 


43.70 


2.98±0.07 


A1795 


0.062 


44.30 


2.92±0.04 


A2029 


0.077 


44.53 


3.07±0.03 


A2065 


0.072 


44.10 


3.03±0.11 


A2142 


0.089 


44.68 


3.05±0.04 


A2256 


0.058 


44.34 


3.13±0.02 


A2319 


0.056 


44.56 


3.19±0.02 


A3112 


0.070 


44.06 


2.74±0.06 


A3266 


0.055 


44.26 


3.04±0.03 


A3376 


0.046 


43.76 


2.84±0.04 


A3391 


0.054 


43.87 


2.82±0.10 


A3395 


0.050 


43.91 


2.93±0.03 


A3558 


0.048 


44.23 


2.99±0.02 


A3571 


0.040 


44.26 


3.02±0.04 


A3667 


0.053 


44.35 


2.99±0.02 


A4059 


0.048 


43.82 


2.93±0.18 


Cygnus A 


0.057 


44.30 


3.20±0.11 


MKW3S 


0.045 


43.73 


2.79±0.04 



a The X-Ray luminosities arc in the 0.1-2.4 keV 
band, from Markevitch (1998). We take uncer- 
tainty in the luminosities to be 20%. 

b From Fadda et al. (1996). 



